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Viscous flows within an elastic structure apply stress on the solid-liquid interface. The stress-field 
created by the viscous flow can be utilized to counter stress created by external forces and thus 
may be applied as a tool for delaying the onset of structural failure. To illustrate this concept we 
study viscous flow within an elastic cylinder under compressive axial force. We obtain a closed- 
form expression showing an approximately linear relation between the critical buckling load and 
the liquid inlet pressure. Our results are validated by numerical computations.We discuss future 
research directions of fluid-solid composite materials which create flow under external stress, yielding 
enhanced resistance to structural failure. 


Viscous flows within an elastic solid apply both pres¬ 
sure and shear stress at the solid-liquid interface and thus 
create an internal stress- and deformation-fields within 
the solid. The stress-field created by the viscous flow 
can be utilized to counter stress-field created by exter¬ 
nal forces, similarly to pre-stressing of solid structures 
[H, and thus may be applied as a tool for delaying the 
onset of structural failure. We illustrate this concept by 
examining the effect of internal viscous flow on buckling 
of an elastic cylinder under compressive axial force. The 
examined case is relevant to applications such as soft- 
robots and soft-actuators UM and micro-needles [3^, 
where the buckling failure mode is one of the limiting 
factors determining the needle radius and thus is of in¬ 
terest for on-going efforts to create less invasive, painless 
micro-needle arrays. 

Various existing works studied the interaction between 
external viscous flow and elastic buckling of a slender fil¬ 
aments. Becker & Shelley (2001) Q studied deformation 
of a slender elastic filament due to a flow-field with uni¬ 
form shear stress and examined bifurcations of the defor¬ 
mation shape. Wandersman et al. (2010) Q conducted 
experiments of elastic fibers within a viscous flow com¬ 
posed of an array of vortices. Guglielmini et al. (2012) Q 
studied, analytically and numerically, the effect of exter¬ 
nal viscous flow on the buckling instability of a filament 
clamped to a rigid wall. Buckling and instabilities of 
elastic cylinders due to internal viscous flow was exten¬ 
sively studied in the context of pipes conveying fluid and 
collapsible tubes [iM3- 

In this work we focus on steady viscous, Newtonian, in¬ 
compressible flow through a slender linearly elastic cylin¬ 
der (see Figure [T]) and study the effects of shear stress 
applied by the liquid on the onset of the buckling failure 
mode due to external axially compressible force. Here¬ 
after, capital letters denote normalized variables and as¬ 
terisk superscript denote characteristic values. We de¬ 
note the length of the cylinder I , the inner radius , the 
outer radius Tq, the Young’s modulus of the solid Cm, the 
critical buckling load fc, the second moment of inertia of 
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the solid cross section the liquid viscosity /i, liquid 
density p, liquid speed in the axial direction w, character¬ 
istic solid radial deformation d*, characteristic solid axial 
deformation d* and characteristic solid stress aT (acting 
on the plane normal to coordinate i and in the direction 
of coordinate j). 

We require small radius-to-length ratio 


I 


= ei < 1 


and small reduced Reynolds number 


„ PW Tj 

Rcr = ei- -<C 1, 
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which is the relevant parameter expressing the ratio of 
inertia to viscosity in slender geometries, where w* is 
the characteristic velocity. We estimate the character¬ 
istic pressure in the liquid by comparing the order-of- 
magnitude of the critical buckling load fe and the total 
axial force applied the by the liquid on the solid 
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p*nrf 
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FIG. 1. A schematic description of the stress applied on the 
cylinder in the unbuckled state (a) and the cylinder cross 
section (b). Part (c) presents the deflection and forces and 
moment acting at 2 ; = Z for the buckled case. 
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where k is effective column length factor and thus p* = 
emmi/pTTrf. 

We define the normalized coordinates 


z r 

Z=-, R=- 

l Vi 


F, = 


fef 


P = 


ppTTrf 


(4) 


corresponding to the normalized axial coordinate, nor¬ 
malized radial coordinate, external force and liquid pres¬ 
sure, respectively. Under the above assumptions the flow 
within the cylinder is Hagen-Poiseuille flow and thus the 
stress applied by the liquid on the solid surface r = is 


AP 

Fizr{Z,R = 1 ) = 2 ) ^rr( 


riZ,R = l) = P{0) + ZAP, 

(5) 

where AP = P(l) — P(0) is the normalized pressure dif¬ 
ference, Yjzr is the dimensionless shear stress applied by 
the liquid, normalized by cr*^ = eip* and Yirr is the di¬ 
mensionless radial stress, normalized by a*^ = p*. Sub¬ 
stituting the Hagen-Poiseuille relation w* = r^p*/8pl 
(where w* is estimated as the average liquid speed) and 
rrii = TTrf{n‘^ — \)/A into ([2]) yields 


« 1, 


8pPl 
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We can thus estimate the characteristic stresses cr*^ 
cr*^, within the cylinder by 


Pirrj 


and 
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’ u — 1 P-Krj{n — 1) 

From the constitutive relation and Hook’s law we obtain 


ddz 

dz 
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and since a*^ ^ a *^, a*r , Ugg we can estimate the defor¬ 
mation field based on as 

rj* o2 rl* rj* 

+ ( 10 ) 

LA Vi I 

and thus the cylinder deformations are negligible for the 
forces and pressures characteristic with buckling. We can 
thus limit our analysis to cases with constant channel 
cross-section. 

We follow the classic approach by Euler (1757) [T^H^ 
and calculate the minimal value of the external force Fg 
for which a solution of the bending equation of the beam 
can be obtained. The governing equation is 
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dz^ 


fe - ApTTrf (^1 - 


w = q(l — z) + m. 

( 11 ) 


We define the normalized parameters 


w* emrriiW* 


,2 73 


M = 


mk‘^1 

emrriiW* 


( 12 ) 


representing, deflection, lateral reaction force at Z = 1 
and moment at Z = 1, respectively. The normalized 
equation is thus 


d^W 

dZ^ 


+ [Pe - AP(1 - Z)]W = Q{1 -Z)-M. (13) 


Substituting C = AP-^/^{Fg - AP(1 - Z)) = 

AP“^/3Pg — Ap3/3(l — Z)) we obtain the Airy equation 
and the homogenous solution Yh is thus 


Yh = CiA(C) + C2P(C) 


(14) 


where A and B are the zero and second Airy func¬ 
tions, respectively. The complete solution can be ob¬ 
tained by variation of parameters as (utilizing the rela¬ 
tion a(c)p'(C)-A'(c)p(C) = iA) 
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The boundary conditions at .Z = 0 are 

„-,z = 0) = 0. = 


( 16 ) 


yielding Ci = C 2 = 0. The boundary conditions at ^ = 1 
can thus be described as functions of Q and M 

W{l)=Tr [ {M + Q-QOx (17) 
Jo 
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These relations enable solving for the minimal value of 
Fg for various boundary conditions and liquid pressure 
differences AP. For the case of fixed end at Z = 1, 
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FIG. 2. The value of k, effective column length factor, as 
function of AP. For all cases the column is fixed a.t Z = Q. 
Ax Z — 1 M = Q and W {1) = 0 for case 1, M = 0 and 
W{1) = 0 for case 2 and dW{\)ldZ = 0 and W(l) = 0 for 
case 3. The smooth lines are solutions of EHli-dlll) and the 
dashed lines are a approximated solution CT- 


dW[X)ldZ = 0 and W{1) = 0, (fT^ and (IT^ are solved 
for Ff, and Q/M simultaneously. For hinged end at Z = 
1, M = 0 and W{\) = 0, (fT51) is solved for Fg. For end 
free to move laterally at Z = 1, Q = 0 and dW{l)/dZ = 
0, dni) is solved for Fe- 

We approximate the obtained solutions to the linear 
expression 

Fe ~ S (1 + CAP) (19) 

/Cq 

or in dimensional form as fc = jkQ^ernTfiijP + 

Cirr^Ap), where fcg is the effective column length fac¬ 
tor calculated for the case without internal flow and C 
is a coefficient representing the effect of the viscous flow. 
The value of C is C = 1/80 for 91F(1)/9Z = 0 and 
W{1) = 0, C = 1/60 for dW(\)ldZ = 0 and M = 0, and 
C = 1/80 for 9FF(1)/9Z = 0 and Q = 0. Equivalently, 
we can present this approximation in terms of effective 
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FIG. 3. An illustration of a simple solid-liquid composite 
material. A small pin with length small compared with the 
length of the cylinder is connected to the base of a liquid- 
filled cylinder and thus flow is created when external force is 
applied. 

lution is evident and a significant monotonic decrease in 
fc, and thus an increase in rigidity, as AP increases. 

Our results can be readily applied as a tool to tem¬ 
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structure which may allow delaying structural failure and 
gaining complex structural deformations. An illustration 
of a simple solid-liquid composite material is presented 
in figure |31 A small pin with length small compared with 
the length of the cylinder is connected to the base of a 
liquid-filled cylinder and thus flow is created when ex¬ 
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